﻿ACADEMIA FORŢELOR TERESTRE “NICOLAE BĂLCESCU” REFERAT PROBLEME DE PROGRAMARE LINIARĂ DE TIP TRANSPORT COORDONATOR Lt col conf univ c s I Vasile CĂRUŢAŞU AUTOR Slt Ioan STERP -SIBIU, 2013- CUPRINS 1 Metode pentru determinarea unei soluţii de bază 1 Metoda minimului pe linie 1 Metoda minimului pe coloană 1 3 Metoda minimului pe tabel 1 4 Metoda colţului de nord-vest 1 5 Metoda diferentelor maxime 2 Metode pentru determinarea unei soluţii optime 2 1 Metoda acoperirii zerourilor 2 2 Metoda circuitelor PROBLEME DE PROGRAMARE LINIARĂ DE TIP TRANSPORT Forma generală a unei probleme de tip transport este dată de tabelul: B1B2 BnDisp A1c11C12 c1nD1 A2c21C22 c2nD2 Amcm1cm2 cmnDm NecesarN1N2 Nm unde A1, A2, , Am sunt firme producătoare care au disponibile cantităţile D1, respectiv D2, , Dm (dintr-un anumit produs), iar B1, B2, , Bn sunt firme care au nevoie de produsul respectiv în anumite cantităţi: N1, N2, , Nm Costul corespunzător transportului unei unităţi de produs de la firma Ai la beneficiarul Bj este dată de costul cij înscris în tabel Dacă totalul disponibilului este egal cu totalul necesarului, adică EMBED Equation 3 , ne punem problema distribuirii cantităţilor de produse disponibile către beneficiari, în condiţiile unor costuri totale de transport minime O soluţie a acestei probleme se scrie sub forma EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 , unde EMBED Equation 3 reprezintă cantitatea de produs pe care o primeşte beneficiarul Bj de la producătorul Ai Pentru a determina costul total de transport, se calculează funcţia obiectiv EMBED Equation 3 1 Metode pentru determinarea unei soluţii de bază Aceste metode duc la găsirea unei soluţii în distribuirea produsului de la producători la beneficiari, fără însă a obţine un cost minim de transport 1 1 Metoda minimului pe linie Se caută costul minim de transport de pe prima linie, EMBED Equation 3 Se compară disponibilul de pe linia 1 cu necesarul de pe coloana j Dacă disponibilul este mai mic decât necesarul, se atribuie celulei x1j disponibilul existent şi se trece la următoarea linie Dacă disponibilul este superior necesarului, atunci se atribuie celulei necesarul corespunzător şi se trece la următoarea celulă cu cost minim (dintre cele rămase) de pe linia 1 Se continuă procedeul până se consumă întregul disponibil de pe linia 1, după care se trece la linia următoare După ce se atribuie întregul disponibil beneficiarilor, se scrie soluţia şi se calculează funcţia obiectiv Exemplu B1B2B3B4Disp A145034550 A21004564100 A3351004504150 A47451006100 Necesar10050100150400 Solutia va fi: EMBED Equation 3 ; functia obiectiv: EMBED Equation 3 1 2 Metoda minimului pe coloană Se caută costul minim de transport de pe prima coloană, EMBED Equation 3 ci1 Se compară disponibilul de pe linia i cu necesarul de pe coloana 1 Dacă disponibilul este mai mic decât necesarul, se atribuie celulei xi1 acest disponibil şi se trece la următorul cost minim de pe coloana 1 Dacă disponibilul este superior necesarului, atunci se atribuie celulei necesarul corespunzător şi se trece la următoarea coloană Se continuă procedeul până se atribuie întregul necesar pentru coloana 1, după care se trece la coloana următoare După ce se atribuie întregul disponibil beneficiarilor, se scrie soluţia şi se calculează funcţia obiectiv Exemplu B1B2B3B4Disp A145034550 A24561004100 A3100355044150 A474505506100 Necesar10050100150400 Solutia va fi: EMBED Equation 3 ; functia obiectiv: EMBED Equation 3 1 3 Metoda minimului pe tabel Se procedează ca la metodele anterioare, cu deosebirea că se caută costul minim din întreaga matrice a costurilor, EMBED Equation 3 Se compară disponibilul cu necesarul şi după încărcarea celulei cu cost minim, se va trece la următoarea celulă cu cost minim din tabel Exemplu B1B2B3B4Disp A145034550 A250456504100 A3503510044150 A47451006100 Necesar10050100150400 Solutia va fi: EMBED Equation 3 ; functia obiectiv: EMBED Equation 3 1 4 Metoda colţului de nord-vest Această metodă nu ţine cont de costul fiecărei celule în alegerea celulei care trebuie încărcată ci de poziţia acesteia în tabel La fiecare pas se va alege celula din colţul de nord-vest (stânga sus) din dreptunghiul rămas Pornim de la celula x11 Dacă disponibilul este mai mare decât necesarul, se atribuie celulei acest necesar şi se trece la celula din colţul de nord-vest a dreptunghiului rămas (acesta este tabelul iniţial fără prima coloană, deoarece întregul necesar pentru beneficiarul B1 a fost atribuit) Dacă disponibilul este mai mic decât necesarul, se atribuie celulei acest disponibil şi se trece la celula din colţul de nord-vest a dreptunghiului rămas (acesta este tabelul iniţial fără prima linie, deoarece producătorul A1 nu mai are nici o unitate de produs disponibilă) Se continuă procedeul până la distribuirea tuturor cantităţilor de produs existente Exemplu B1B2B3B4Disp A150434550 A250450564100 A3351004504150 A47451006100 Necesar10050100150400 Solutia va fi: EMBED Equation 3 ; functia obiectiv: EMBED Equation 3 1 5 Metoda diferentelor maxime Pe fiecare linie şi pe fiecare coloană a tabelului se calculează diferenţa dintre cel mai mic cost de transport şi cel imediat superior Se identifică linia sau coloana cu cea mai mare diferenţă şi aici, în celula de cost minim, se execută prima alocare Se refac diferenţele pe liniile şi coloanele neblocate folosindu-se numai costuri "neblocate", după care se reia procedura de alocare Exemplu B1B2B3B4Disp A145034550 A24561004100 A3100354504150 A47410056100 Necesar10050100150400 Solutia va fi: EMBED Equation 3 ; functia obiectiv: EMBED Equation 3 2 Metode pentru determinarea unei soluţii optime Aceste metode conduc la o soluţie optimă din punct de vedere al costului transportului Cele mai cunoscute metode de determinare a soluţiei optime sunt metoda acoperirii zerourilor şi metoda circuitelor(ciclurilor) 2 1 Metoda acoperirii zerourilor Într-o primă etapă se urmăreşte obţinerea a cel puţin unui zero pe fiecare linie şi coloană din matricea costurilor În acest scop , se scade elementul minim al fiecărei linii din toate elementele liniei respective şi, după epuizarea tuturor liniilor, se scade elementul minim al fiecărei coloane din toate elementele coloanei respective În a doua etapă se construiesc toate acoperirile posibile O acoperire este alcătuită din două mulţimi I şi J; I este o mulţime de linii ale matricei costurilor (cu maxim m-1 linii), iar J este o mulţime de coloane, formată din acele coloane care conţin zerouri, mai puţin pe liniile din I Astfel se consideră prima linie şi toate coloanele ce conţin zerouri în afara primei linii Se consideră pe rând fiecare linie cu coloanele ce conţin zerouri în afara liniilor considerate, iar apoi se construiesc acoperirile cu 2,3, ,m-1 linii Pentru fiecare acoperire (I,J) astfel construită se efectuează suma dintre cantităţile disponibile aferente liniilor din I şi cantităţile necesare corespunzătoare coloanelor din J Dacă toate sumele calculate sunt mai mari sau cel puţin egale cu echilibrul problemei , se trece la etapa a treia În caz contrar , se consideră acoperirea cu suma cea mai mică şi în matricea costurilor se taie cu o linie orizontală fiecare linie a acoperirii şi cu linie verticală fiecare coloană a acoperirii Astfel se obţin trei tipuri de elemente: elemente netăiate, elemente tăiate o singură dată şi elemente tăiate de două ori Se consideră elementul minim netăiat şi el se scade din toate elementele netăiate Totodată acest element minim netăiat este adăugat la elementele dublu tăiate Elementele tăiate o singură dată rămân neschimbate Pentru noua matrice astfel obţinută se construiesc din nou toate acoperirile posibile şi se reia astfel etapa a doua până când toate sumele calculate devin mai mari sau cel puţin egale cu echilibrul problemei , moment în care se trece la etapa a treia În cazul în care în noua matrice nu s-a pierdut nici unul dintre zerourile existente anterior, atunci este suficient să se studieze doar acele acoperiri cărora le corespunseseră sume mai mici decât echilibrul problemei, deoarece , nepierzând nici un zero, sumele calculate nu pot scădea, ci doar eventual să crească În etapa a treia se încarcă soluţia în poziţiile în care se află zerouri, începând cu acea linie sau coloană care conţine cel mai mic număr de zerouri Se continuă în acest mod încărcarea până când se epuizează tot disponibilul, iar toată cererea din punctele de consum a fost satisfăcută Exemplu B1B2B3B4Disp A1434550 A24564100 A33544150 A47456100 Necesar10050100150400 Matricea costurilor este: EMBED Equation 3 Scadem elemental minim de pe fiecare linie si apoi pe coloane Obtinem: EMBED Equation 3 Trecem acum la etapa de acoperire a liniilor Pentru EMBED Equation 3 se observă că zerourile din afara liniei 1, se găsesc pe coloanele 1,2,3,4, deci EMBED Equation 3 Suma EMBED Equation 3 o calculăm adunând disponibilul existent pe linia 1 cu necesarul existent pe coloanele 1, 2, 3, 4, adică EMBED Equation 3 Procedăm analog pentru toate celelalte combinaţii de linii Obţinem EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 Trecem acum la combinaţiile de câte două linii: EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 Cum combinaţiile de două linii au fost epuizate trecem acum la acoperirea a câte trei linii: EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 Dacă toate sumele obţinute pentru fiecare din acoperirile scrise mai sus ar fi mai mari sau cel mult egale cu totalul disponibilului (sau al necesarului) - în cazul exemplului considerat – 400, atunci s-ar trece la încărcarea soluţiei optime În caz contrar, se caută toate sumele mai mici decât 400 şi se alege dintre acestea cea cu valoarea cea mai mică Am subliniat mai sus toate variantele cu S<400 Avem un singur astfel de caz: EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 Modificăm acum matricea costurilor în felul următor Taiem liniile si coloanele corespunzatoare: Obtinem noua matrice: EMBED Equation 3 In noua matrice nu mai exista zeroul de pe pozitia (1,2), deci vom reface toate acoperirile Recalculand sumele pentru noua matrice se va observa ca toate sunt mai mari sau egale cu 400, deci am obtinut optimul Vom trece la incarcarea solutiei Zerourile din matricea costurilor modificată indică celulele din tabelul iniţial care pot fi încărcate cu diverse cantităţi pentru a obţine soluţia optimă căutata Matricea este EMBED Equation 3 , deci celulele ce pot fi incarcate sunt: EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 , EMBED Equation 3 Zerourile din matricea costurilor modificată indică celulele din tabelul iniţial care pot fi încărcate cu diverse cantităţi pentru a obţine soluţia optimă căutata Procedand asa cum este prezentat in descrierea algoritmului, incarcarea tabelului va fi: B1B2B3B4Disp A145034550 A24561004100 A3100354504150 A47410056100 Necesar10050100150400 Solutia va fi: EMBED Equation 3 ; functia obiectiv: EMBED Equation 3 2 2 Metoda circuitelor Într-o primă etapă se determină o soluţie iniţială nedegenerată (o soluţie care să aibă un număr de componente nenule egal cu numărul producătorilor plus numărul beneficiarilor minus 1) utilizând una din metodele cunoscute În etapa a doua, pentru fiecare celulă (i,j) rămasă neocupată , deci pentru care EMBED Equation 3 , se construieşte un circuit care să aibă toate vârfurile doar în căsuţe încărcate, deci pentru care EMBED Equation 3 Pentru fiecare circuit astfel construit se calculează o valoare EMBED Equation 3 astfel: începând cu căsuţa neocupată de la care s-a început construirea ciclului respectiv, se adună şi se scad succesiv costurile unitare de transport corespunzătoare colţurilor circuitelor Dacă toate valorile astfel obţinute sunt mai mari sau egale cu zero, atunci soluţia este optimă şi algoritmul se încheie Dacă însă există şi valori negative, soluţia poate fi îmbunătăţită Astfel, se alege circuitul cu cea mai mică valoare şi în colţul liber al său se va încărca o cantitate fictivă t, iar în celelalte colţuri , cantitatea existentă va fi modificată alternativ cu ( t Se va determina t drept cantitatea maximă cu care se poate modifica soluţia astfel încât componentele ei să rămână nenegative În acest mod a fost determinată o nouă soluţie căreia i se aplică din nou algoritmul descris la etapa a doua; algoritmul se încheie atunci când toate valorile ataşate circuitelor sunt mai mari sau egale cu zero Daca nici una din metodele pentru determinarea unei solutii de baza nu duce la obtinerea unei solutii nedegenerate, sau daca degenerarea apare pe parcursul aplicarii metodei atunci, in momentul aparitiei degenerarii se va adauga la una din liniile sau coloanele care duc la degenerare o cantitate infima z Pentru pastrarea echilibrului problemei, aceeasi cantitate infima z va fi adaugata si la o linie/coloana care inca nu a fost epuizata Se aplica algoritmul ca mai sus, iar la final, dupa obtinerea optimului, se aloca lui z valoarea 0 Se va obtine astfel solutia optima a problemei initiale B1B2B3B4Disp A1434550 A24564100 A33544150 A47456100 Necesar10050100150400 Exemplu Incepem prin a gasi o solutie de baza nedegenerata Am obtinut deja solutiile de baza pentru aceasta problema si toate cele 5 metode au dus la obtinerea de solutii degenerate In aceasta situatie vom alege una dintre metode Alegem metoda coltului de nord-vest Solutia obtinuta a fost: B1B2B3B4Disp A150434550 A250450564100 A3351004504150 A47451006100 Necesar10050100150400 Degenerarea a aparut la incarcarea celulei (2,2), cand s-au epuizat si linia a doua si coloana a doua Vom adauga o cantitate infima z fie la linia 2 fie la coloana 2 Sa adunam la linia 2 Pentru pastrarea echilibrului vom adauga z si la una dintre coloane Nu vom adauga la coloana 2, deoarece degenerarea s-ar pastra Vom adauga la coloana 3 Noul tabel, cu solutia oferita de metoda coltului de nord-vest va fi: B1B2B3B4Disp A150434550 A2504505Z64100+z A3351004504150 A47451006100 Necesar10050100+z150400 Pentru fiecare din celulele necompletate din tabel se formează câte un circuit format numai din celule completate, schimbând semnul la fiecare schimbare de direcţie a circuitului Pentru celula (1,2) avem circuitul: EMBED Equation 3 În continuare: EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 Pentru fiecare cij de mai sus, se calculează costul total EMBED Equation 3 al circuitului, scriind în relaţiile de forma EMBED Equation 3 , costul existent în fiecare celulă a circuitului (semnele rămânând aceleaşi) EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 EMBED Equation 3 Dacă toate valorile calculate EMBED Equation 3 ar fi mai mari sau egale cu zero, soluţia din tabel ar fi cea optimă În cazul problemei considerate, există mai multe valori negative pentru EMBED Equation 3 , ceea ce înseamnă că soluţia trebuie modificată Pentru aceasta ne vom opri la cea mai mică dintre valorile negative EMBED Equation 3 În cazul problemei, valoarea minimă este -2, corespunzătoare mai multor circuite Alegem circuitul: EMBED Equation 3 In colţul liber al său se va încărca o cantitate fictivă t, iar în celelalte colţuri , cantitatea existentă va fi modificată alternativ cu ( t Se va determina t drept cantitatea maximă cu care se poate modifica soluţia astfel încât componentele ei să rămână nenegative B1B2B3B4Disp A150-t43T4550 A250+t4505z-t64100+z A3351004504150 A47451006100 Necesar10050100+z150400 Vom avea EMBED Equation 3 (retinem ca z este o cantitate infima, mai mica decat toate celelalte) Noua solutie astfel obtinuta va fi: B1B2B3B4Disp A150-z43z4550 A250+z450564100+z A3351004504150 A47451006100 Necesar10050100+z150400 Se refac circuitele pentru noua solutie si se verifica optimul Repetand de cateva ori algoritmul, vom obtine valori nenegative pentru toate costurile atunci cand solutia va fi: B1B2B3B4Disp A143504550 A2456100+z4100+z A310035Z450-z4150 A475045056100 Necesar10050100+z150400 Vom aloca lui z valoarea 0, solutia optima a problemei fiind data de tabelul: B1B2B3B4Disp A143504550 A24561004100 A3100354504150 A475045056100 Necesar10050100150400 Avem: EMBED Equation 3 si EMBED Equation 3 